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SECTION  I 

INTRODUCTION  AND  PROJECT  GOALS 


This  report  describes  the  goals,  the  progress  to  date  in  two  areas  of 
effort,  and  the  plans  for  future  work  on  Contract  NObsr  91332.  The  four 
major  goals  on  this  project  are  listed  below.  The  first  two  are  of  primary 
importance  since  the  latter  two  cannot  be  fully  achieved  without  reaching  the 
first  two. 

^ The  development  of  refined  methods  of  “Automated  LOFAR-gram 
reading"  is  an  improvement  to  the  PACT  techniques  of  extracting  spectral 
clues  from  a received  signal.  It  is  now  intended  that  the  implementation  of 
these  techniques  be  primarily  by  a general-purpose  digital  computer.  This 
approach  will  allow  greater  freedom  in  the  choice  of  processing  techniques, 
will  avoid  the  problems  of  interconnecting  the  existing  equipment  to  a new 
type  of  computer,  and  will  permit  a more  thorough  investigation  to  be  car- 
ried out  within  the  available  funds. 

Development  of  practical  methods  of  predicting  (computing)  impor- 
tant spectral  characteristics  of  passive  sonar  signals  from  a knowledge  of 
the  generating  ship's  (target's)  engineering  and  operating  configuration. 

More  specifically,  this  problem  is  approached  by  developing  generic 
models  of  certain  vibratory  ship's  components.  Integral  to  the  task  of  de- 
veloping prediction  methods  are  the  subtasks  of  collecting  the  two  types  of 
data  required  — tape  recordings  of  ship-generated  acoustic  signals  and  engi- 
neering data  on  the  corresponding  ships. 

Development  of  a classification  algorithm  which  operates  on  the 
spectral  clues  extracted  by  the  LOFAR-gram  reader,  obtaining  a set  of 
probable  ship  design  characteristics,  and  thence  a set  of  probable  target 
classifications.  

Implement  the  procedures  developed  above  in  a manner  suitable  for 
use  in  further  studies  and  for  test  and  evaluation  of  these  techniques. 
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SEC  TION  II 

SIGNAL  PROCESSING  STUDIES 


This  section  describes  a study  of  spectrum  analysis  in  passive  sonar 
classification,  including  the  results  obtained  to  date. 

2.  1 BACKGROUND  AND  REASONS  FOR  THE  STUDY 

Passive  sonar  signal  processing  has  traditionally  made  extensive  use 
of  spectrum  analysis  for  detection  and  classification  purposes.  Spectrum 
analysis  techniques  will  be  used  extensively  by  Automatic  Passive  Sonar  Tar- 
get Classification  (APSTAC)  and  will  be  combined  in  a later  study  with  transi- 
ent processing  to  provide  a complete  classification  methodology. 

Current  plans  call  for  APSTAC  signal  processing  to  be  performed  on 
a general-purpose  digital  computer  (using  analog -to -digital  conversion  and 
digital  tape  storage  of  the  signals).  One  of  the  key  factors  making  this  approach 
possible  was  the  recent  development  of  the  Fast  Fourier  Transform  (FFT)^ 
which  dramatically  reduces  the  time  and  thus  the  cost  of  computer  spectrum 
analysis.  The  FFT  is  a computational  algorithm  which  reduces  the  number 
of  time  consuming  operations  required  in  transforming  a set  of  equispaced  data 
points . 

The  study  described  in  this  progress  report  was  performed  to  generate  a 
definite  base  for  the  analysis  of  real  passive  sonar  signal  recordings,  and  con- 
sisted of: 

a.  Generation  of  the  required  computer  programs 

b.  Checkout  of  these  programs 

c.  Familiarization  with  the  properties  of  these  programs  through 
simple  examples 

d.  Investigation  of  the  methodologies  of  computer  spectrum  analysis 
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e.  Investigation  of  proposed  methods  of  signal  processing  using  com- 
puter simulations  of  passive  sonar  or  derived  signals 

f.  Investigation  of  the  effects  of  quantization  and  sampling  rate  on 
the  effectiveness  of  the  signal  processing 

NOTE;  This  study  provided  an  opportunity  for  gaining  further 
insight  into  the  influence  of  the  FFT  on  general  computing 
methods. 

2.2  COMPUTER  METHODS  OF  SPECTRUM  ANALYSIS 

Digital  computers  may  be  used  to  perform  spectrum  analyses  of  time 
waveforms  by  each  of  several  different  approaches.  Until  recently,  the  ac- 
cepted method  was  to  compute  the  appropriate  auto-  or  c ross -correlation  func- 
tion (mean-lagged  product)  for  delays  less  than  some  fraction  of  the  total  ana- 
lyzed data  length,  and  then  to  compute  a fraction  of  the  Fourier  coefficients  of 
the  correlation  function.  This  method  yielded  only  a small  fraction  of  the  avail- 
able spectral  information,  and  was  thus  time  consuming  and  costly. 

A series  of  recent  papers^  ’ ^has  described  a new  algorithm  for  com- 
puting Fourier  coefficients  which  is  dramatically  faster,  for  practical-sized 
problems,  than  the  algorithms  used  previously.  As  a result  the  most  efficient 
method  of  computing  spectral  data  is  now  to  Fourier  transform  the  entire  data 
sequence  to  be  analyzed,  compute  the  squared  magnitudes  of  the  Fourier  co- 
efficients, and  perform  smoothing,  if  appropriate.  In  addition,  the  FFT  also 

(4) 

now  provides  the  fastest  method  to  compute  correlations  or  convolutions. 

Except  for  overhead  functions,  the  time  required  by  the  FFT  to  compute  a set 

of  T Fourier  coefficients  is  proportional  to  T log  T whereas  the  computation 

2 

times  for  the  algorithms  used  previously  were  proportional  to  T . 

A digital  computer  operates  on  signals  which  have  been  quantized  in 
both  amplitude  and  time.  An  independent  theoretical  study  is  now  under  way 
at  Litton  to  determine  the  effect  on  the  degree  of  accuracy  in  amplitude  quan- 
tization with  respect  to  the  fidelity  of  the  spectra  obtained.  Results  of  this 
other  study  will  be  incorporated  into  APSTAC  as  soon  as  they  become  available. 
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Figure  2-1.  Frequency  Axis  Folding  to  Produce  Aliasing 
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The  effect  of  time  sampling  on  the  computed  spectrum  may  be  des- 
cribed as  shown  in  Figure  2-1  and  is  termed  aliasing.  The  power  which  is 
at  the  frequency  nfg  ±f,  where  f = 1/A  is  the  sampling  rate  and  n is  an  in- 
teger, cannot  be  distinguished  from  power  at  frequency  f.  Thus,  the  trans- 
formation of  the  continuous  frequency  spectrum  to  the  discrete  frequency  spec- 
trum may  be  thought  of  as  a folding  operation  as  shown  in  Figure  2-1  followed 
by  sampling  in  frequency.  (It  should  be  noted  that  fg/2  is  often  called  the  fold- 
ing frequency.  ) Thus,  the  sampling  error  in  estimating  spectra  depends  upon 
the  rate  at  which  the  spectrum  falls  off  with  increasing  frequency.  If,  for  ex- 
ample, the  highest  frequency  of  interest,  B,  is  f /4  and  the  spectrum  falls  off 
at  a rate  of  6dB  per  octave,  then  the  aliasing  power  will  be  about  9dB  below 
the  desired  power  level. 


The  FFT  operates  on  a finite  sequence  of  complex  data,  x(t)  for  t = 0, 
1,  2,  ...,  T-l,  and  yields  the  finite,  discrete  Fourier  transform,  X (y)’ 
of  the  time  sequence  as: 


T-l 


*(t) 


Y J x(t)e'lZ,r  kt/T’  for  k = 0,  1,  2, 


. . , T-l 


t=0 


The  inverse  transform  of  X (y)  *s; 


T-l 


x(t)  = 


I 

k=0 


k \ i2ff  kt/T 


* (t)* 


Throughout  most  of  this  section  the  time  variable  t has  been  normalized  to 
an  integer.  Real  time  will  then  be  tA,  where  A is  the  sampling  interval. 
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which  is  continuous,  periodic  with  period  T,  and  takes  on  the  values  $?(t)  = x(t) 
for  t = 0,  1,  2,  ....  T-l.  This  equality  may  be  demonstrated  as  follows.  For 
t = 0,  1 T-l, 


x(t)  - Sc(t)  = x(t)  - J X (±) 


k \ i2ff  kt/T 

rr-  1 ^ 


5 *(t>  - J 


T-l  T-l 

v i2ir  kt/T  1 V / » -i2ir  kv/T 


T S x(v)  e 


= x(t,  - 2 x(v)l  £ 


= x(t)  - x(t)  = 0 


where  the  orthogonality  condition: 


4i2ir(m  - k)t/T 


T if  k = m 


0 if  k * m 


has  been  used. 

To  show  that  the  complex  exponential  is  orthogonal  over  a finite  set  of 
equispaced  points,  note  that: 


el2ir  tl/T  = ^ ei2ff  U/T  for  1 s < s T - 1 
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so  that  the  sum  of  the  first  forward  differences  of: 


/ti\  _ iZit  t-f /T 

It)  - e 


is  zero.  That  is, 


2 -s-®  [•(*)-• 

t=0  1 t=0 


i2lM/T  , 
= e - X 


T-l 

2 


= 0 for  1 s /.  s x - 1 


Thus,  the  summed  product  of  complex  exponentials  is  zero  unless  the  argu- 
ments are  the  same.  That  is, 


2 


2ir  a/T  =0  if  l * 0 (mod  T) 


The  other  theorems  which  are  conventionally  stated  for  continuous  time 
and  frequency  variables  also  carry  over  to  the  discrete,  finite  sample  point 
case.  For  example,  Parseval's  theorem  holds;  the  sum  of  the  powers  at  each 
frequency  is  equal  to  the  total  power  in  the  time  series.  That  is. 
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= I -4  2 *,„  * w2' k<<  - t)/T 


k=0  * t,  1=0 


- i I *<•>  2. 


where  an  asterisk  is  used  to  denote  complex  conjugation. 

The  finite,  discrete  Fourier  coefficients  have  several  important  prop- 
erties. The  coefficients  are  periodic  with  period  T. 


j x(t)e-i2,,,<ktT»/T  = X(t) 


T T 

Therefore,  the  coefficients  for  k = y+  1,  tj-  + 2,  ....  T-l  are  those  that 

are  normally  associated  with  negative  frequencies. 

(k  \ sjc 

Yj  for  simplicity  of  notation.  If  x(t)  is  real,  a^  = a 

and  so, 


ak  + a.k  = ak  + aT.k  = 2Re  ak 


Because  of  the  periodicity  of  the  sequence  a^  , if  x(t)  is  real  and  T is  even, 


x(t)  = £ 


i2tr  kT/T 


- a + /a  eiZ*t/T  + a ei2l7t/T\  + + a 

" ao  + \al  C + aT-l  e j + ...  + aT/2 
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2fft  , 
cos  — + 


I a I 
I a2l 


cos 


4irt 


* 2Re  laT/2-l|  cos  2l,UIfZ 


^ * <-‘>‘aT/2 


Similarly,  if  x(t)  is  imaginary  and  T is  even, 


T/2-  1 

x(t)  = aQ  + 2 £ Im  ]a^|  cos + (-l)4  aT/2 

-t=l 


Convolution  in  the  time  domain  corresponds  to  the  dual  operation  in  the 
frequency  domain  if  the  convolution  is  performed  "circularly;"  i.  e.  , as  if  the 
time  sequences  were  periodic: 


T-l 

^ x(t)  y(t  - T ) = y 
t=0 


Now  that  some  of  the  more  important  properties  of  the  finite  discrete 
Fourier  transform  have  been  summarized,  it  remains  to  be  decided  as  to  how 
the  Fourier  coefficients  are  to  be  used  in  computing  a "power  spectral  density" 
or  "spectrum.  " The  theory  and  practice  of  spectrum  computation  and/or  esti- 
mation are  undergoing  continual  development.  The  general  subject  is  one  of 
considerable  complexity  because  of  the  many  facets  which  must  be  considered. 
The  method  to  be  used  in  any  particular  case  should  be  chosen  with  regard  to 
the  specific  conditions  and  objectives  of  that  case. 

If  one  is  computing  a spectrum  of  a time  function  which  is  characterized 
by  a stable,  strong  line  structure  it  is  usually  important  to  minimize  the  "leak- 
age" of  power  at  a frequency  f * k/T  for  any  integer  k;  i.  e.  , f is  not  any  of 
the  sampled  frequencies,  over  to  the  spectral  estimates  at  frequencies  removed 
from  f.  The  most  widely  accepted  way  of  reducing  leakage  is  to  "hann"  the 
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data,  either  by  multiplying  the  time  series  by  C 1 - cos  2irt/T]  (a  cosine  bell), 
or  by  convolving  the  Fourier  coefficients  with  the  three  element  sequence  -1/2, 

1,  -1/2.  Hanning  the  data  by  either  method  changes  the  power  leakage  from  a 
| f — k/ T | ^ dependence  co  a |f  - k/T|  ^ law.  ^ 

If  the  time  waveform  to  be  spectrum  analyzed  is  random  with  a fairly 
flat  spectrum,  one  is  faced  with  the  necessity  of  trading  off  some  frequency 
resolution  for  increased  statistical  stability  of  the  estimated  spectrum.  This 
can  be  achieved  by  locally  smoothing  the  Fourier  coefficients  by  convolving 
with  a suitably  chosen  short  sequence. 

The  simplest  type  of  spectrum  which  can  be  computed  from  the  Fourier 
coefficients,  X(k/T),  is  the  (unmodified)  periodogram,  P(k/T),  given  by: 

P (k/T)  = |X  (k/T)  |2 

Tnis  is  a natural  definition  for  a spectrum,  but  one  which  has  drawbacks  in 
many  cases.  The  operations  of  hanning  and  smoothing  of  the  Fourier  coef- 
ficients as  described  above  result  in  modified  periodograms.  For  the  purposes 
of  this  study  to  date,  however,  the  spectrum  computed  has  been  the  unmodified 
periodogram. 

2.  3 SIMULATION  OF  PASSIVE  SONAR  SIGNALS 

During  the  course  of  the  APSTAC  study,  a number  of  different  methods, 
involving  spectrum  analysis,  of  processing  passive  sonar  signals  will  be  studied 
both  theoretically  and  experimentally.  In  order  to  test  these  methods  of  signal 
processing  under  completely  known  conditions  with  controlled  degrees  of  vari- 
ability, several  programs  have  been  written  which  generate  time  series  simu- 
lating the  properties  of  passive  sonar  signals  which  are  currently  of  interest; 
i.  e.  , harmonic  groups  in  a noise  background.  The  programs  have  been  written 
in  a general  manner  so  that  input  parameters  such  as  the  number  of  harmonic 
groups,  the  fundamental  frequencies,  the  amount  of  phase  drift  of  each  har- 
monic group,  the  amplitudes  of  the  individual  harmonics,  etc.  , can  be  speci- 
fied with  ease. 
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In  addition  a simple  noise  generation  routine  has  been  prepared  which 
can  serve  as  the  background  noise  for  the  harmonic  groups,  or  as  a signal 
source  for  studying  the  effects  of  sampling  rate.  The  routine  written  is  basi- 
cally a one-pole,  low-pass  filtering  operation  on  the  output  of  a normal  or 
gaussian  random  number  generator  (white  noise).  In  continuous  time,  such 
a filter  has  an  impulse  response: 


h(t)  = v/ia  e'at  U^t) 


where  U^t)  is  a unit  step  at  t = 0 and  a is  a parameter  which  specifies  the 
filter  effective  smoothing  time, 


H(f) 


+ 


and  has  an  asymptotic  falloff  of  6dB  per  octave. 

The  digital  equivalents  of  the  above  continuous  functions  are: 


h^t)  = y 1 - e"2a  e'at  6 1(t)  U^t) 

and 

Hi<« 

where  &a(t)  is  a periodic  unit  impulse  train  with  period  a,  and  * denotes  con- 
volution. This  form  of  the  transfer  function  indicates  that  the  one-pole  char- 

- 1 

acteristic,  (1  + i 2trf/a)  , is  repeated  at  unit  intervals  along  the  frequency 
axis,  and  thus  graphically  indicates  the  restrictions  on  a to  avoid  aliasing. 

The  half-power  bandwidth  of  the  one -pole  filter  is  Bj  = a/2ir  and  for 
an  nth  order  cascade  is  Bn  = Bj  y/zl^n  - i.  Therefore,  n = 2,  B2  = 0.  643B1 
and  for  n = 3,  B^  = 0.  51Bj. 
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The  actual  implementation  of  this  filter  on  a computer  is  recursive  in 
form.  That  is,  if  r(t),  t = 0,  1,  ...  is  the  input  sequence,  a new  sequence  is 
formed  by 


P(t)  = r(t)  + e'ap(t-l) 

and  this  is  then  scaled  by  A - e 2a  to  retain  unity  power  transfer. 

If  the  input  to  this  filter  r(t),  is  a sequence  of  zero  mean,  independent 

2 

random  numbers  (white  noise),  with  variance  a then  the  autocovariance  func- 
tion of  the  output,  r(t),  is 

C?  ?(t,  t-T)  = J|?(t)?(t-T)| 


E 


r(|)e'(t'i)a  a]>  r(m)e 


-(t-T  -m)a 


1 = 0 


m=0 


= ar2e-aT  [l  - (e-2a)t-T  + 1] 


The  bracketed  term  displays  the  effect  of  the  starting  transient  on  the  filter 
output.  After  the  filter  "fills  up,  " however,  the  output  may  be  said  to  possess 
a correlation  function  (dependent  only  upon  the  shift  variable  t),  i.e.  , 


Ry,  y(t) 


lim  C-  -(t,  t-T)  = a 2 e"aT 

t-a,  r*  r r 


Hence,  the  filter  must  be  initialized  by  several  iterations  before  the  output 
is  added  as  background  to  the  simulated  harmonic  group  structures. 

The  first  signal  generation  subroutine  to  be  described  here,  SIGGEN  2 
(S1GGEN  2 is  a generalization  of  an  earlier  subroutine  called  SIGGEN  1),  is 
intended  to  simulate  a direct  passive  sonar  signal  by  a sum  of  harmonic  groups 
plus  a background  noise.  Up  to  three  harmonic  groups  can  be  generated  and 
each  group  has  a specific  set  of  parameters: 
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Fj’  Aj’  nj’  Ajl’  Aj2’  • 


where  n.  s 11.  The  form  of  the  jth  harmonic  group  is 


where  9-n  is  a uniform  [0,  2ir]  random  number,  independent  over  both  j and  n. 
T values  of  sj(t)  are  generated  in  a block  as  an  input  to  the  spectrum  analysis 
routine.  Each  harmonic  group  has  its  own  independent  random  phase  drift, 
Cpj(t)  which  is  generated  by 


9j«> 


0 if  t = 0 

\j0j(t)  + cph(t-l)  for  t > 0 


where  i/).(t)  is  a sequence  of  random  numbers  independent  over  both  t and  j 
J 1/2 

and  uniform  [_aj>  ajl>  with  a^  = (3FjT)  . This  choice  of  a^  results  in  an 
rms  phase  drift  of  radians  per  cycle.  The  output  of  this  routine  can,  at 
option,  be  quantized  into  a specified  number  of  levels.  Examples  of  spectra 
of  signals  generated  by  this  routine  will  be  given  in  the  next  subsection. 


A subroutine,  S1GGEN  3,  has  been  written  which  produces  an  output 
similar  to  that  observed  when  a passive  sonar  signal  is  spectrum  analyzed. 
The  spectrum  is  read  out  serially,  in  time,  and  the  spectrum  as  a time  func- 
tion is  one -bit  quantized  on  a dynamic  threshold  which  is  just  above  the  con- 
tinuous portion  of  the  spectrum.  This  series  of  operations  is  illustrated  in 
Figure  2-2.  Such  a binary  signal  may  be  used  for  detecting  the  harmonic 
groups  and  other  line  structure  properties  of  a passive  sonar  signal.  (Just 
such  a signal  was  used  as  the  input  to  the  PACT  correlator.  ) 


The  output  of  S1GGEN  3 is  a binary  sequence  which  is  a sum  of  signal 
terms  simulating  the  effect  of  harmonic  groups  in  the  hypothetical  input  pas- 
sive sonar  signal,  and  a noise  term  which  simulates  the  effect  of  noise  spikes 


2-11 

CONFIDENTIAL 


CONFIDENTIAL 


i 


in  the  continuous  portion  of  the  "input  spectrum"  exceeding  the  one -bit  quanti- 
zation threshold.  The  probability  of  an  "input  noise  spike"  exceeding  the  one- 
bit  quantization  threshold  is  assumed  to  be  constant  over  the  entire  system 
bandwidth  so  that  the  probability  of  a noise  pulse  in  the  simulation  output  is 


n(t)  = 


0 with  probability  1 - g 

1 with  probability  (3 


for  t = 0,  1,  2 T-  1.  g is  an  input  parameter  to  the  program. 

The  signal  components  in  the  SIGGEN  3 output  consist  of  roughly  peri- 
odic. binary  pulse  trains,  where  the  probability  of  pulse  occurrence  decays 
exponentially  with  increasing  time  (time  corresponds  to  frequency  in  the  sonar 
spectrum).  Let  the  interval 
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The  form  of  an  individual  signal  component  is  then 


0 if  trfl. 

s.(t)  = ^ 0 with  probability  1 - P.e  if  tfL^ 

1 with  probability  P.e  ajm  if  tfL^ 


for  t = 0,  1,  . . . , T - 1,  and  m = 1,  2,  . . . [(T  - 1 )/ (D^  + d./2)]. 

There  can  be  up  to  three  signal  terms  and  for  each  there  is  an  input  parameter 
array 


D.,  d.,  P.,  a. 
J J J J 
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The  output  binary  sequence  is  generated  from  the  signal  and  noise 
terms  by  the  summing  operation 


!0  if  y(t)  = ^ s (t)  + n(t)  = 0 
j=l 

1 if  y(t)  a 1 

where  ISR  is  the  number  of  signal  terms.  Examples  of  spectra  of  signals 
generated  by  SIGGEN  3 will  be  given  in  the  next  subsection. 

The  SIGGEN  subroutines  are  combined  with  a fast  Fourier  transform 
subroutine,  called  FFT,  and  an  output  subroutine,  called  OUTPP,  to  make 
up  the  functional  portions  of  a general  software  package  called  the  Frequency 
Domain  Processing  Program  (FDPP).  The  general  operations  of  the  FDPP 
are  diagrammed  in  Figure  2-3.  This  program  has  been  written  in  FORTRAN 
and  used  on  an  IBM  360/50  computer. 

2.  4 EXAMPLES  OF  SPECTRA  OF  SIMULATED  SIGNALS 

In  order  to  illustrate  the  properties  of  the  finite  discrete  Fourier  trans 
form,  a special  computer  run  was  made  using  SIGGEN  1.  The  time  function 
which  was  spectrum  analyzed  was 


-i2ir  e~r?  t 80 

512  + cos  2 if  g j-j 


80  , 150.  5 

5TT  t + cos  — 5TZ“  t 


. 30°-  25  * 4.  9-  511.5  . 

+ cos  2tr — — 1 + cos  t 

for  t = 0,  1,  2,  . . . , T = 512,  which  includes  a complex  term  and  four  real 
terms.  The  complex  time  function  has  a spectral  component  only  at  f^  = -70/512 
whereas  a sine  or  cosine  of  the  same  frequency  would  have  spectral  components 
at  ±fj.  Two  of  the  terms  in  §(t)  are  at  frequencies  of  the  form  f = n/T  where 
n is  an  integer.  These  frequencies  are  among  those  sampled  by  the  finite, 
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Figure  2-3.  Frequency  Domain  Processing  Program  1734-13 

Flow  Chart 


discrete  Fourier  transform  and  these  terms  do  not  affect  the  values  of  the 
other  spectral  components  (with  the  exception  of  computational  roundoff). 
(This  is  due  to  the  orthogonality  of  the  complex  exponentials,  or  sines  and 
cosines,  over  a finite,  equispaced  set  of  points.  ) 

In  contrast  to  the  effect  of  terms  at  frequencies  of  the  form  f = n/T 
where  n is  an  integer,  terms  at  frequencies  of  the  form  f = X/T  where  X is 
not  an  integer,  affect  the  values  of  every  spectral  component.  Naturally, 
terms  in  §(t)  at  frequencies  of  the  form  f = X/T  contribute  most  strongly  near 
their  own  frequencies  and  their  effect  on  other  spectral  components  falls  off 
rapidly  with  |k  - X | / T in  a manner  dependent  upon  how  close  X is  to  an  inte- 
ger. All  these  effects  are  shown  in  Figure  2-4  which  is  a plot  of  the  periodo- 

jJc 

gram  of  ?(t). 

sjc  sjc 

5(t)  contains  two  terms  at  frequencies  above  the  Nyquist  rate,  namely 
at  the  frequencies  300.  25/ 512  and  511.  5/ 512.  A term  in  5(t)  above  the  Nyquist 
frequency  shows  up  at  (or  around)  its  low  frequency  aliase,  in  this  case 
1 - (300.  25/ 5 12)  = (2 1 1 . 75/ 512)  and  (0.  5/ 512),  respectively.  Therefore,  the 


All  the  spectra  presented  here  have  been  simplified  to  an  extent.  The  spectra 
are  functions  of  a discrete  variable  k,  but  portions  of  them  have  been  plotted 
as  if  they  were  continuous.  Furthermore,  the  spectra  to  follow  have  been 
smoothed  (visually)  over  the  low  power  density  portions. 

Since  time  has  been  normalized  to  the  integers,  the  Nyquist  rate  here 
is  1/2. 
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Figure  2-4.  Spectrum  of  Five  (Complex)  Sine  Waves 


1734-14 


spectral  components  shown  at  k = ±211.  75  and  at  k = ±0.  5 in  Figure  2-4  show 
the  effect  of  undersampling,  i.e.,  sampling  at  less  than  twice  the  frequency 
of  terms  in  the  time  function  which  have  significant  power.  Note  that  the  zero 
frequency  term  is  not  strongly  aliased  by  cos  2ir(511.  5/512)t,  rather,  the  value 
shown  at  k = 0 is  primarily  due  to  the  term  in  §(t)  with  frequency  150.  5/512. 
Computational  roundoff  error  is  negligible  in  all  of  the  described  runs. 

As  described  in  the  previous  subsection,  the  subroutine  S1GGEN  2 
generates  signals  which  are  harmonic  groups  with  random  phase  drifts.  Fig- 
ure 2-5  shows  a (smoothed)  spectrum  of  such  a signal.  The  analyzed  time 
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Figure  2-5.  Spectrum  of  Three  Harmonic  Groups  with  1734-15 

Random  Phase  Shift 

function  consists  of  the  sum  of  three  signal  groups,  each  of  which  consists  of 
a cosine  at  the  fundamental  frequency  of  the  group  plus  a cosine  at  the  seventh 
harmonic  of  the  group  fundamental.  There  is  no  noise  background  to  these 
harmonic  groups. 

The  first  signal  term,  at  the  fundamental  frequency  10.  5/1024,  has  no 
random  phase  drift  and  so  its  spectral  line  widths  are  essentially  the  same  as 
similar  lines  in  Figure  2-4.  The  seventh  harmonic  is  84/1024,  a sampled 
frequency,  and  so  the  harmonic  line  is  narrower  than  the  fundamental  spec- 
tral line. 
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The  second  signal  term,  at  a fundamental  frequency  of  30.  5/  1024,  has 
an  rms  phase  drift  of  0.  104  radians  per  cycle.  The  width  of  the  fundamental 
line  is  about  the  same  as  would  be  the  case  without  phase  drift,  but  the  width 
of  the  seventh  harmonic  line  at  244/1024  is  substantially  greater  when  phase 
drift  is  present. 

The  third  signal  term,  at  a fundamental  frequency  of  50.  5/1024,  has 
an  rms  phase  drift  of  0.  190  radians  per  cycle.  The  width  of  the  spectral  line 
at  the  fundamental  frequency  is  unaffected  by  the  random  phase  drift.  The 
seventh  harmonic,  however,  is  greatly  spread  in  frequency  and  the  local  peak 
of  the  power  density  does  not  occur  at  8(50.  5/1024)  = 404/1024,  but  rather  the 
local  peak  power  is  at  408/1024. 

As  shown  in  Figure  2-5,  it  can  be  concluded  that  in  terms  of  simulating 
passive  sonar  signals,  the  rms  phase  drift  per  cycle  of  the  fundamental  must 
be  maintained  less  than  one -tenth  radian. 

The  signal  generation  subroutine  SIGGEN  3 outputs  a binary  sequence 
which  is  similar  to  the  signal  used  for  harmonic  group  detection  in  the  PACT 
feasibility  hardware.  As  an  initial  illustration  of  the  type  of  spectra  possessed  by 
signals  generated  by  SIGGEN  3,  Figure  2-6  is  the  periodogram  of  a nonrandom 
train  of  slowly  widening  pulses.  (Relatively  little  smoothing  has  been  performed 
on  the  spectrum  in  Figure  2-6.  ) More  precisely,  the  time  series  which  re- 
sulted in  the  spectrum  in  Figure  2-6  may  be  represented  by  a pulse  train  with 
Period  7,  the  first  40  pulses  being  one  unit  wide,  and  the  remaining  33  pulses 
being  3 units  wide.  No  pulses  are  missing  and  there  are  no  extra  (noise ) pulses . 
In  terms  of  the  signal  parameters  defined  in  Subsection  2.  3,  the  signal  here  is 
determined  by 


T = 512;  D - 7;  d = 0.  05;  P = 1;  a = 0 

Therefore,  the  spectrum  of  this  signal  is  a set  of  harmonically  related  lines 
with  fundamental  frequency  T/D  = 73,  as  seen  in  Figure  2-6. 
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of  Pulse  Occurrence 


Although  no  conclusions  can  be  drawn  from  Figure  2-7  regarding  the 
efficacy  of  spectral  analysis  as  a technique  for  detecting  harmonic  groups,  at 
least  a partial  illustration  of  the  limitations  of  such  techniques  is  given. 
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SECTION  III 

CLASSIFICATION  DATA/SHIP  MODELING  PROGRESS 

3.  1 FUNCTIONAL  DESCRIPTION  OF  A FIRST  VIBRATION/ RADIATION 

MODEL 

Modeling  of  a main  propulsion  system  as  the  principal  exciter  of  the 
so-called  line  spectrum  of  acoustic  radiations  is  well  underway  under  the 
present  contract  effort.  Objectives  of  the  initial  effort  have  been  to  identify 
or  develop,  where  necessary,  mathematical  descriptors  of  the  harmonic  con- 
tent of  energy  inputs,  frequency  response,  and  energy  transmission  charac- 
teristics of  a particular  ship  propulsion  system. 

Several  considerations  influenced  the  selection  of  a ship  for  model 
development.  The  selected  ship  should  be  a relatively  uncomplicated  one  for 
which  detailed  engineering  data  is  available.  Further,  it  is  desirable  that  all 
data  be  unclassified  to  simplify  collection  and  handling.  The  basic  propulsion 
system  should  give  rise  to  a model  which  can  be  generalized  to  fit  naval  ships 
of  high  interest. 

These  considerations  led  to  the  modeling  of  a single-screw,  diesel- 
driven  merchant  ship  of  a design  submitted  to  the  U.S  Maritime  Administra- 
tion by  the  J.J.  Henry  Company,  naval  architects,  of  New  York.  The  model 
of  this  ship  can  readily  be  modified  to  permit  sound  field  syntheses  of  many 
classes  of  merchant  ships  as  well  as  diesel -driven  naval  ships  such  as  sub- 
marines, certain  classes  of  ships;  i.  e.  , mine-sweepers,  etc. 

The  approach  taken  in  the  vibrational/acoustic  radiation  modeling  of 
ships  is  a conservative  one.  That  is,  ships  are  finite  systems  and  consider- 
able knowledge  exists  regarding  the  vibrational-radiative  mechanisms  by  which 
their  radiated  sound  fields  are  produced.  The  synthesis  of  composite  models 
will,  of  course,  rely  heavily  on  this  knowledge  fund,  some  of  which  is  based 
on  substantial  investigation  and  years  of  effort. 

The  engineering  vibrational  characteristics  of  ship  systems  in  which 
we  are  most  interested  are  those  which  cannot  be  eliminated  nor  reduced 
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Figure  3-1.  Ship's  Acoustic  Vibration/Radiative  System 
Configuration 
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below  some  level  which  is  a function  of  the  ship's  operating  power  at  any 
given  time.  The  main  propulsion  system,  including  the  engines,  gearing, 
shafting  and  propellers,  is  the  principal  source  of  this  type  of  excitations. 
Certain  auxiliary  machines  also  contribute  excitations  which  can  prove  use- 
ful but  these  are  not  incorporated  into  the  present  model  since  their  acoustic 
evidences  are  usually  separable  from  those  of  the  propulsion  train. 

3.  1.  1 Elements  of  the  Problem 

Evolution  of  a first  systematic  model  of  a ship  propulsion  system  in- 
volved several  problems  which  are  best  visualized  with  reference  to  the  dia- 
gram shown  in  Figure  3-1.  Each  function  included  in  the  propulsion  system 
line  spectrum  synthesizer  under  construction  either  contributes  a prolific 
spectrum  of  input  energies  or  serves  as  a complex  frequency-sensitive  filter 
of  those  excitations.  The  problems  of  developing  a composite  model  are 
principally: 
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a.  To  evolve  mathematical  descriptors  which  are  practically 
manipulable. 

b.  To  generate  outputs  from  each  function  which  are  compatible 
inputs  to  those  which  follow. 

c.  To  account  for  all  dispositions  of  energy  from  input  to  output 
and  to  retain  all  phase  relationships  through  the  sequence  of 
calculations. 

d.  To  develop  plausible  descriptors  for  vibrational/radiative  mecha- 
nisms which  are  empirically  known  to  be  present  and  functional 
but  for  which  mathematical  descriptors  are  not  available  in 
literature. 

e.  To  devise  algorithms  on  which  to  base  calculations  for  computer - 
programming  and  data  manipulations  in  complex  format,  to  or- 
ganize the  various  routines  into  a system,  to  make  provision  for 
examining  results  at  several  points  within  the  model  and  to  pre- 
sent line  spectra  and  other  outputs  in  their  most  useful  form. 

The  succeeding  paragraphs  discuss  the  nature  of  these  problems  and 
the  means  employed  to  date  for  their  solution. 

3.  1.  2 General  Description  of  the  Basic  Model 

The  basic  model  is  a computer  program  which  synthesizes  the  vibra- 
tional excitation  of  a marine  diesel  propulsion  system.  This  excitation  divides 
its  energy  between  two  radiators;  the  propeller  and  the  ship's  hull.  The  pro- 
peller is  simultaneously  an  influential  element  of  the  excitation  system  and  a 
radiator.  It  is  natural  to  expect  that  these  two  roles  are  interdependent  and 
that  any  meaningful  mathematical  description  of  the  propeller  as  an  inertia 
and  damper  in  the  excitation  system  should  be  affected  by  its  role  as  an  acoustic 
projector.  Thus,  a model  of  a ship's  propeller  is  needed  which  reconciles  the 
two  functions.  The  basic  model  has  evolved  from  considerations  such  as  these, 
and  is  derived  from  engineering  practice  in  the  vibration  analysis  of  marine 
propulsion  systems  and  from  the  need  to  merge  diverse  analytical  procedures 
into  a unified  computer  program. 
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The  general  nature  of  the  model  which  has  been  developed  can  be  de- 
scribed from  one  standpoint  by  stating  the  principle  assumptions  on  which  it 
is  based.  First,  it  is  assumed  that  the  propulsion  system,  as  an  exciter,  can 
be  modeled  as  an  integral  system  of  forcing  functions,  passive  elements  and 
propeller -radiator  so  long  as  radiation  reaction  effects  at  the  propeller  are 
accounted  for  and  that  similar  effects  from  the  hull  do  not  significantly  in- 
fluence the  characteristics  of  the  exciting  system.  This  does  not  preclude 
the  hull  excitation  nor  radiation.  This  assumption  states  that  it  is  not  an  in- 
fluential element  of  the  excitation  system.  Figure  3-2  illustrates  the  system 
treated  in  the  present  model. 

It  is  further  assumed  that  the  propulsions  system  can  be  modeled  as 
a linear,  lumped  parameter,  mechanical  system.  Linearity  in  this  case 
means  that  the  vibrational  results  of  forcing  functions  applied  at  the  individual 
cylinders  will  combine  by  superposition  at  the  output.  Lumping  of  parameters 
is  not  performed  where  available  engineering  data  permits  them  to  be  described 
separately  and  where  doing  so  will  make  the  model  more  definitive  and  infor- 
mative. The  reciprocating  components  and  the  rotating  components  of  each 
cylinder  assembly  are  separated  and  their  effects  calculated.  However,  the 
compliances  of  coupled  sections  of  shafting  between  inertias  are  lumped  be- 
cause it  does  not  degrade  the  ability  of  the  model  to  serve  its  purpose. 

The  mobility  electromechanical  analogy  is  taken  as  a useful  technique 
for  representing  this  propulsion  system.  It  results  in  an  electrical  network 
of  the  form  shown  in  Figure  3-3  representing  a 9-mass,  8-compliance  system 
with  three  included  dampers.  Dampers  are,  of  course,  the  only  devices  by 
which  energy  can  be  converted  or  lost  in  a mechanical  system.  The  most 
important  damper  in  this  system  is  the  propeller  which  converts  a steady 
torque  (analogous  to  direct  current)  to  the  thrust  which  drives  the  ship  and 
converts  torsional  vibrations  (analogous  to  ac  electrical  power)  to  waves 
of  excess  pressure;  i.  e.  , to  acoustic  radiations.  It  should  be  noted  that  the 
network  representation  in  Figure  3-3  is  principally  for  convenience  in  setting 
up  equations  in  terms  which  will  be  immediately  meaningful  to  a majority  of 
engineers. 
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Figure  3-3.  Analog  Representation  of  a Marine  Diesel  1734-4 

Propulsion  System 

The  pressure  inputs  which  occur  within  the  cylinders  are  obtained 
from  actual  indicator  diagrams  made  of  the  cylinder  pressures  of  a marine 
diesel  engine.  Such  pressure  curves  are  rich  in  harmonics,  and  the  form  for 
inputting  the  data  from  the  curves,  as  sample  points,  is  sufficient  to  permit 
harmonic  analysis  up  to  some  desired  number  of  harmonics.  Calculations 
are  carried  out  to  convert  the  sample  values  from  input  pressures  to  input 
torques  along  the  crankshaft.  This  is  both  economical  of  computation  effort 
and  generates  a harmonic  torque  series  which,  in  the  mobility  analogy,  is 
represented  by  electrical  current.  Thus,  the  network  presents  complex  im- 
pedances to  harmonic  series  of  input  currents  at  as  many  points  as  there  are 
cylinders  in  the  engine.  The  input  series  for  the  various  cylinders  differ  only 
by  a phase  angle  determined  by  the  number  of  pistons  and  the  firing  order. 

It  should  be  noted  that  the  selection  of  firing  order  is  an  important  design  con- 
sideration in  the  vibration  control  of  marine  diesel  propulsion  systems. 
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The  solution  of  the  network  for  values  of  the  input  impedance  it  pre- 
sents at  each  cylinder  for  all  harmonics  of  input  current  (torque)  is  compli- 
cated by  the  fact  that  the  magnitudes  of  damping  in  the  system  are  functions 
of  frequency.  The  impedance  values  must  be  computed  separately  for  each 
running  speed  of  the  engine  and,  in  the  case  of  propeller  damping,  for  each 
frequency  of  harmonic  output.  The  initial  form  of  output  of  the  network  is  a 
series  of  voltages  across  the  elements  representing  the  propeller.  That  is, 
superposition  of  effects  is  assumed  and  phasor  summation  of  the  outputs  due 
to  each  current  input  is  made,  on  a harmonic -by-harmonic  basis.  Voltage 
in  the  mobility  analogy  represents  velocity,  a series  of  sinusoidal,  torsional 

vibration  velocities  in  the  present  case.  To  calculate  power  outputs  the  ana- 

y2 

log  relationship  P = is  used,  recalling  that  both  V and  R have  separate 
values  for  each  harmonic  of  output. 

The  development  of  a simple  mathematical  model  of  a propeller  as  an 
acoustic  projector,  which  is  consistent  with  the  hydrodynamic  model  of  a pro- 
peller, and  which  provides  reasonably  accurate  predictions  of  signature  char- 
acteristics, will  be  a significant  and  far-reaching  contribution.  It  is  believed 
that  the  results  to  date  (through  preliminary)  have  been  significant  and  prom- 
ising. The  specific  assumptions  on  which  this  part  of  the  overall  model  is 
based  are  stated  in  the  discussions  of  Subsection  3.  1.  3. 


3.1.3  Functional  Components  of  the  Basic  Model 

Figure  3-4  is  a functional  flow  diagram  of  the  marine  diesel  propulsion 
system  model  now  available  for  investigations  of  the  vibrational-accoustical 
radiations  of  a specific  class  of  ships.  This  diagram  is  also  the  gross  flow 
chart  of  the  computer-programming  which  implements  the  model.  To  review 
the  functions  shown  in  Figure  3-4: 

a.  Blocks  I and  IX  read  in  data  points  taken  from  two  harmonically 
rich  functions.  The  cylinder  pressure  indicator  curve  for  a particular  engine 
is  assumed  common  to  all  the  cylinders  of  that  engine.  Subsequent  products 
formed  on  this  set  of  data  points  support  an  harmonic  (Fourier)  analysis,  the 
results  of  which  are  an  harmonic  set  of  torque  inputs  to  the  engine.  The  pro- 
peller wake-inflow  data  used  are  curves  for  (1  - Wj),  the  weighted  (from  root 
to  tip)  mean  longitudinal  pressure  wake  encountered  through  2n  radians  of 
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rotation  by  a propeller  blade.  (A  typical  curve  is  subsequently  shown  in 
Figure  3-6  in  this  report.  ) As  many  as  243  data  points  from  each  curve  may 
be  supplied  the  Fourier  analysis  routine  permitting  up  to  121  harmonics  of 
content  to.be  extracted. 

b.  Block  II  converts  the  indicator  pressure  data  points  for  a cylinder 
to  a corresponding  set  of  torque  values  at  the  crank.  That  is: 

2 

Qk  = Pk  R sin  (^),  k = 1,2 K 

where  K is  the  number  of  pressure  data  points,  D is  the  piston  diameter, 
are  the  pressure  values,  and  R is  the  crank  throw  of  the  engine. 

c.  Blocks  III  and  X employ  the  same  finite  Fourier  analysis  routine 

to  generate  the  Fourier  coefficients,  a and  b , for  harmonic  series  of  N 
6 n n 

terms,  the  number  of  harmonics  of  propeller  shaft  rotation  frequency  for 
which  the  program  will  predict  vibrational  and  acoustical  outputs. 

d.  Due  to  the  varying  force  arm  in  the  torque  applied  by  a piston  in 
a reciprocating  engine  as  the  crank  rotates  through  2n  radians  of  arc,  an 
harmonically-varying  inertial  torque  is  generated  which  must  be  added  to  the 
harmonic  driving  torque.  Though  actually  an  infinite  harmonic  series  itself, 
the  inertial  torque  is  adequately  accounted  for  by  phasor  addition  of  its  first 
three  terms  to  the  series  computed  in  Block  III  so  long  as  the  L/R  ratio  (see 
below)  is  large  enough.^  This  is  accomplished  in  Block  IV.  The  terms  added 
are  described  by: 

Q.  ...  = 1/  2M  oj  ^R^  ( sin  ut  - sin  2tut  - sin  3 cut  ) 

inertial  rec  \ 2L  2L  / 

Where  Mrec  is  the  mass  of  the  reciprocating  components  of  one  piston  assem- 
bly, o>  is  the  radian  frequency  equal  to  the  propeller  shaft  rotation  frequency, 
R is  engine  crank  throw  as  before,  and  L is  the  length  of  the  connecting  rod, 
usually  4 to  5 times  R. 


3-9 

CONFIDENTIAL 


CONFIDENTIAL 


e.  Blocks  V and  XI  add  phase  angles  to  the  harmonic  torque  series. 

Block  V adds  phase  to  account  for  the  lag  due  to  firing  order  of  the  input  torques 
at  each  cylinder.  Block  XI  accounts  for  the  phase  displacement  of  the  propeller 
blades  from  a reference  blade  in  rotating  through  the  variable  wake  inflow.  It 
also  adds  an  arbitrary  phase  angle  ^between  Je  ) accoun*  for 

the  fixed  difference  in  angular  position  of  the  reference  propeller  blade  and  the 
number  one  engine  crank. 

f . Block  VII  operates  on  the  network  representation  of  the  propulsion 
system  (see  Figures  3-2  and  3-3).  Input  impedances  are  computed  for  all 
harmonics  of  input  current  (for  reasons  given  in  the  next  paragraph)  at  each 
of  the  input  nodes  (cylinder  locations).  Output  voltages  at  all  harmonic  fre- 
quencies across  the  output  resistance  (propeller  damping)  are  obtained  by 
phasor  summation  of  contributions  due  to  the  inputs  from  each  cylinder. 

The  necessity  for  computing  separate  network  impedances,  both 
input  and  output,  for  all  harmonic  frequencies  arises  from  the  fact  that  the 
propeller  functions  as  an  acoustic  projector.  This  is  intrinsic  to  its  function 
as  a torque-to-thrust  converter  and  to  pertinent  hydrodynamic  considerations 
which  are  well  developed  in  the  literature  related  to  propeller  design  and  the 
powering  of  ships.  ^ A unified  description  of  the  propeller  as,  simultaneously, 
a torque -to -thrust  converter,  as  an  acoustic  projector,  and  as  a variable  (with 
frequency)  damper  and  energy  storage  device  has  been  developed  for  use  in  the 
present  model. 

g.  Blocks  VI  and  VII  compute  two  contributions  which  combine  to 
make  up  the  all-important  propeller  damping  and  a variability  (with  respect  to 
frequency)  of  the  polar  mass  moment  of  inertia  of  the  propeller.  The  viscous 
damping  of  propeller  torsional  vibration  at  a frequency  equal  to  the  propeller 
shaft  rotation  rate,  which  differs  from  the  propulsive  resistance  to  the  steady- 
state  propeller  rotation,  can  be  obtained  in  two  ways.  One,  is  from  the  usual 
rpm  versus  speed  and  shaft  horsepower  versus  speed  curves  of  the  ship. 
Assuming  power  P,  torque  Q,  and  shaft  speed  0 given  in  basic  units,  they  are 
related  by: 

P = Q6  or  Q = P/d 
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Damping,  Kq,  is  given  by  dQ/d 0 so  that: 

K = i dP/d0  - P/02 

° 0 

- i.  _P_  pounds /foot 

0 \dS  / \dl?/  qZ  radians  / second 

where  S is  ship's  speed  in  knots. 

Now  dP/dS  and  dS/d0  are  just  slopes  which  can  be  taken  from  the 
curves  mentioned  above  for  the  appropriate  ship's  speed. 

The  second  way,  and  one  which  permits  accounting  for  variables 
such  as  propeller  slip  assumed  constant  in  the  preceding  calculation,  is 
through  use  of  an  extensive  set  of  curves  worked  up  by  Archer^  based  on 
model  experiments  and  empirical  data.  These  curves  yield  a value  for  the 
coefficient  a to  be  used  in  the  equation 


K = a Q /6 
o 


pounds / foot 
radians  /second 


The  two  methods  yield  closely  similar  values  for  Kq.  Block  VI  is  an  attempt 
to  tabulate  or  otherwise  store  in  the  computer  Archer's  data  in  such  form  that 
interpolations  into  the  data  for  variations  in  propeller  slip,  pitch-diameter 
ratio  and  disc -area  ratio  can  be  programmed  as  part  of  the  overall  model. 

h.  The  torsional  vibration  damping  discussed  in  the  last  paragraph 
is  applicable  only  to  vibration  at  the  fundamental  frequency.  It  is  well  known 
that  the  damping  is  greater  for  higher  harmonic  frequencies.  Den  Hartog^^ 
suggests  simply  doubling  the  value  of  damping  obtained  by  computing  dQ/d# 
previously  described. 

The  more  scientific  approach  is  to  seek  a model  of  the  propeller 
as  an  acoustic  projector  and  to  observe  the  results  of  radiation  reaction  on 
the  propeller.  The  radiation  reaction  on  a projector  is  a complex  quantity, 
the  real  part  of  which  adds  to  the  damping  and  the  imaginary  part  of  which 
effectively  adds  to  the  mass  moment  of  inertia  of  the  projector. 
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Recalling  that  as  a torque -to -thrust  converter,  a propeller  im- 
parts momentum  to  a column  of  water  astern.  The  reaction  to  this  action  is 
forward  momentum  acquired  by  the  ship.  Thus,  an  ideal  propeller  will  pro- 
ject water  only  in  a direction  directly  opposite  to  tha*  in  which  it  is  desired 
that  the  ship  move.  Actually,  the  column  of  water  is  projected  somewhat 
downward  for  hydrodynamic  reasons,  but  propeller  designs  in  general  do  at- 
tempt to  achieve  the  cylindrical  column  of  wake  immediately  behind  the 
propeller . 

However,  practical  design  considerations  degrade  the  performance 
of  a propeller.  There  must  be  a hub  at  the  center,  and  the  blades,  near  their 
roots,  must  be  designed  more  for  strength  than  for  their  efficiency  in  torque- 
to-thrust  conversion.  Accepting  the  outer  two-thirds,  radially,  of  a propeller 
blade  as  the  principal  region  over  which  torsional  vibrations  at  the  hub  are 
converted  to  acoustic  pressure  trains  projected  aft,  a model  emerges  of  the 
outer  blade  which  is  approximately  equivalent  to  a vibrating  circular  piston. 

This  is  the  key  assumption  on  which  calculations  of  added  damping 
and,  incidentally,  of  added  polar  mass  moment  of  inertia  at  the  propeller  in 
the  present  model  are  based.  The  propeller  is  then  represented  as  an  array 
of  circular  pistons  as  illustrated  in  Figure  3-5.  The  radius  of  an  equivalent 
vibrating  piston  (EVP)  is  expected  to  be  about  one -third  the  route -to -tip 
length  of  a blade  with  the  EVP  tangential  to  the  tip  circle  of  the  propeller. 

Experimental  data  from  actual  ships  will  be  needed,  of  course,  to 
establish  precisely  what  radius  of  circular  piston  will  represent  a blade  of  an 
actual  propeller  or  whether  a more  complicated  shape  would  improve  the  model, 
such  as  perhaps  an  ellipse.  Block  VII  computes  total  damping  at  the  propeller 
for  each  harmonic  frequency, 


K = 
n 


K + K 


rn 


Kq  was  discussed  in  the  preceding  paragraph  and  is  obtained  by  integrating 
for  radiation  damping  over  an  EVP  as  described  above. 
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rn 


poCBR1  (2kna)P 

Tn 


rra 


(r  - a)  + 


Si 


Where 


p o is  the  mass  density  of  the  medium  (sea  water  in  the  present  case) 
C is  the  speed  of  sound  in  the  medium 
B is  the  number  of  propeller  blades 
P is  the  pitch  of  the  propeller 
r is  the  radius  of  the  propeller 

a is  the  radius  of  the  equivalent  vibrating  piston 

k = uu  . is  the  wave  number 
n n/  c 

uun  is  the  radian  frequency  of  the  nth  harmonic,  and 


R.  (2k  a) 
1 n 


(2kna)2 

2-4 


(2k  a)4  (2k  a)6 

n , n 

7 7 

2-4.6  2*4^*  6*8 
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Block  VII  also  computes  total  polar  mass  moment  of  inertia  (PMMI)  of  the 
propeller  including  added  PMMI  due  to  radiation  reactance  at  the  propeller 
blades  and  also  added  PMMI  due  to  entrained  water.  The  expressions  for 
these  quantities  are: 


J = J + J + J 
n o rn  w 


Jq  is  the  PMMI  due  only  to  propeller  mass  and  geometry. 


Jr„  = F X1  <2V>  j 


na.  » 2 . .2 

— «—  - n a (r  - a) 
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where 


X (2k  a)  = l 
in  tt 


2k  a (2k  a) J (2k  a) 
n n 

32-  5 32  • 52-  7 


Jw  = (0.  25)  (PDR)  Jq 


where  PDR  is  the  propeller  pitch/diameter  ratio. 

The  series  Rj  (2k^a)  and  (2k^a)  are  the  piston  resistance  func- 
tion and  piston  reactance  function  originally  derived  by  Lord  Rayleigh  and  de- 
scribed in  his  "Theory  of  Sound.  " 

i . Block  XII  combines  the  torsional  vibration  excitation  transmitted 
via  the  propeller  shaft  to  all  blades  with  the  harmonic  wake  inflow  pressure 
variations  experienced  by  each  blade  as  it  rotates  through  2 n radians  of  arc. 
This  becomes  a polynomial  multiplication  of  the  torsional  vibration  series, 
analogous  to  output  voltages  in  our  network  representation,  with  the  harmonic 
wake  series  which,  it  will  be  shown,  is  equivalent  to  harmonic  variation  of 
the  value  of  the  network  output  resistance. 

It  is  evident  that  harmonic  analysis  of  a mean  longitudinal  wake 
curve,  as  illustrated  in  Figure  3-6,  will  yield  a constant  term  which  is  the 
average  value  of  the  curve  over  the  2n  radians  for  which  it  is  plotted.  This 
could  be  expressed  as 


(1  -w4) 


where 


w-)  = (1  - w#)  + ^ (sinusoidal  components  up  to  N harmonics). 
1 1 n=l 
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Figure  3-6.  Mean  Longitudinal  Wake  Inflow  and  Tangential 
Wake  Curves 
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Repeating  the  equations  for  propeller  damping  at  these  harmonic  frequencies 
we  have 


K = K + K 
n o rn 


poCBR,  (2k  a)  P (_  2 


The  factor  P,  propeller  pitch,  as  normally  used,  is  a constant- 
valued statement  of  a design  characteristic  which  is  the  advance  to  be  ex- 
pected in  one  revolution  of  the  propeller  in  the  absence  of  any  slip.  However, 
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there  always  is  slip;  i.  e.  , failure  to  advance  the  design  distance  per  revo- 
lution. The  generally  accepted  view  is  that  due  to  the  mass  of  a ship,  the 
slip  in  response  to  propeller  torsional  vibrations  is  a function  of  frequency 
but  not  of  time  under  steady-state  operating  conditions,  and  its  value  ap- 
proaches 100  percent  at  high  frequencies.  Slip,  S,  pitch,  and  wake  inflow 
are  related  by  the  expression 


P 


Vs0  - w.) 

-W-TTT 


where  Vg  is  the  ship's  speed  and  6 is  the  shaft  rotation  state. 

Thus,  the  factor  P as  used  in  the  expression  for  K should  not 

rn 

be  thought  of  as  the  constant  but  rather  the  ratio  on  the  right  whose  numerator 
includes  the  wake  inflow  variability  with  time  and  the  denominator  of  which 
includes  the  variability  with  frequency  due  to  slip.  It  is  especially  significant 

to  regard  slip  in  the  above  expression  as  an  element  of  the  factor  j-. =-^  . As 

slip  approaches  100  percent  this  factor  should  become  asymptotic  to  a very 
large  number.  In  order  to  use  a constant  value  for  P in  the  expression  for 

radiation  damping  (assuming  constant  V and#),  the  variabilities  of  (1  - W, ) 

1 s 1 
and  ^ must  be  separately  accounted  for.  Upon  closer  examination,  it 

appears  that  the  quantity  (2k^a)  describes  a variability  with  frequency  quite 
similar  to  that  expected  of  1/(1  - S)  and  for  much  the  same  physical  reasons. 
Plots  of  Rj  (2k^a)  and  (Zk^a),  normalized  by  the  variable  2k^a,  are  pre- 
sented in  Figure  3-7.  This  leaves  only  the  variability  of  (1  - W|)  unexpressed 

in  K as  it  stands.  This  variability  is  described  by  the  harmonic  series  which 
rn 

resulted  from  Fourier  analysis  of  the  wake  inflow  curve,  which,  in  effect, 
results  in  varying  time  and  frequency  damping.  The  reason  for  frequency 
series  multiplication  in  Block  XII  becomes  clear.  Damping,  in  the  mobility 
analogy,  is  the  inverse  resistance  or  conductance.  The  analog  voltage  (vi- 
brational velocity)  across  the  conductance  from  Block  VIII  has  already  been 
obtained.  Output  current  is  now  the  product  of  voltage  and  conductance.  Having 
voltage  and  current,  output  power  is  readily  computed. 
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Figure  3-7.  Piston  Radiation  Reaction  Functions 


1734-2 


1 

I 

1 

0 

0 


It  should  be  noted  that  the  frequency  series  multiplication  is  of  two  com- 
plex polynomials  which  in  the  present  model  may  have  up  to  121  terms  each. 
Their  product  will  consist  of  2N-1  terms,  or  as  many  as  241.  Each  product 
term,  in  turn,  resolves,  using  trigonometric  identities,  into  two  terms  one 
at  the  sum  and  one  at  the  difference  of  the  two  frequencies  entering  into  the 
particular  product  term.  Only  sum  and  difference  terms  at  frequencies  not 
higher  than  those  of  the  two  input  series  are  preserved.  That  is,  only  the  first 
121  of  the  241  terms  resulting  from  the  multiplication  are  regarded  as  useful 
because  significant  contributions  to  terms  122  through  241  could  be  expected 
which  are  not  generated  by  this  multiplication. 

j.  The  mathematics  and  computer  algorithms  for  Block  XII  have  not 
yet  been  developed.  For  the  present,  a simple  rms  summation  of  the  output 
power  computed  for  each  blade  of  the  propeller  is  performed.  This  power  is 
expressed  in  watts  for  each  harmonic  frequency  up  to  the  121st  harmonic  of  a 
fundamental  at  the  frequency  of  propeller  shaft  rotation. 
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k.  Block  XIV,  now  being  set  up,  will  yield  an  output  magnetic  tape 
which  can  be  taken  to  a CALCOMP  plotter.  This  will  in  turn  generate  plots 
of  output  power  spectra  resembling  those  of  Figures  3-8,  3-9  and  3-10.  The 
vertical  scale  will  be  acoustic  output  power  in  dB  relative  to  10”^  watts. 

3.  2 EXERCISING  OF  THE  BASIC  MODEL 

The  two  purposes  of  exercising  the  first  model  are  to  gain  insights  into 
the  vibrational  behavior  of  the  ship  it  represents  and  to  perfect  the  model  as 
an  investigative  tool.  As  stated  earlier,  the  basic  model  which  is  presently 
being  used  is  of  the  central  excitation  system  which  leads  a ship  to  radiate  its 
so-called  line  spectrum  of  propulsion-associated  frequencies.  Whereas  the 
effort  to  date  synthesizes  the  principal  excitation,  work  remains  to  be  accom- 
plished in  describing  the  composite  sound  field.  It  is  important,  however,  that 
the  excitation  mechanism  be  thoroughly  understood  as  the  hull  array  and  the 
combining  with  the  propeller  and  hull  radiations  in  the  near  field  are  described. 

The  term  "exercising  the  model"  is  used  to  represent  the  procedure  of 
varying  parameters  within  the  program  mathematics  to  observe  the  effects  par- 
ticular parameters  have  on  the  output  power  spectral  density.  In  the  case  of 
parameters  which  can  be  varied  independently,  exercising  consists  of  stepping 
those  segments  of  the  computer  program  in  which  the  parameter  operates 
through  solutions  of  their  functions  through  some  range  of  parameter  values. 

To  date  the  model  has  been  exercised  by  introducing  data  points  for  three  dif- 
ferent wake-inflow  curves  into  the  program.  These  are  the  (1  - w^)  curves 
in  Figure  3-6  for:  (1)  the  Mariner  stern,  Propellers  II  and  VI,  (2)  the  conven- 
tional stern,  Propellers  II  and  VI,  and  (3)  conventional  stern,  Propeller  I.  The 
power  spectral  output  differences  corresponding  to  these  differences  in  the 
basic  ship  are  seen  in  the  plots  of  Figures  3-8,  3-9  and  3-10,  respectively. 

It  is  evident  that  exercises,  such  as  those  described  above,  involve 
large-scale  substitutions  of  data  and  may  call  for  special-purpose  computer 
routines  to  analyze  the  effects.  A routine  is  being  written  which  will  subtract 
a trend  from  the  whole  set  of  spectral  outputs  and  perform  a harmonic  analysis 
on  the  residuals  to  ascertain  the  extent  of  harmonic  groupings  within  the  output 
spectrum. 
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Figure  3-9.  Peak  Power  Spectrum  — Output  2 1734-9 
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Exercises  of  interest  to  be  carried  out  using  the  present  basic  model 
may  be  grouped  into  two  categories,  and  at  least  two  analyses  requiring  their 
own  computer  routines  can  be  listed. 


Exercises  to  examine  effects  of  varying  parameters  peculiar  to 
diesel  drives 


(1)  Changes  of  cylinder  firing  order 


(2)  Changes  in  distribution  of  damping  over  the  length  of  the  drive 
train 


(3)  Changes  of  ship's  speed  within  a range  normal  for  the  ship 

(4)  Variation  of  polar  mass  moments  of  inertia  of  fly  wheel,  in- 
tegral and  gear-driven  auxiliaries 

(5)  Attachment  of  devices  such  as  a Bibby  detuner  to  the  engine 

(6)  Variation  of  fixed  phase  difference  between  engine  crank  and 


propeller  blades 


(7)  Variation  of  engine  crank-throw  to  length  of  connecting  rod 
ratio,  R/L 


Exercises  to  examine  effects  of  varying  parameters  common  to 
propellers  and  ships  in  general 


(1)  Differences  in  configuration  of  stern  structure 

(2)  Differences  in  propeller  design,  including  number  of  blades 

(3)  Differences  in  ship  loading  and  resistance  due  to  wind,  sea 
state,  etc.  , as  reflected  in  change  of  propeller  slip 

(4)  Variations,  within  ranges  consistent  with  resistance-powering 
considerations,  of  propeller  pitch -diameter  ratio  and  disc- 


area  ratio 


(5)  Variation  of  assumptions  relating  an  equivalent  vibrating  pis- 
ton to  actual  configuration  of  a propeller  blade 
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c.  Analyses  utilizing  computer  routines 

(1)  Holzer  method  analysis  to  determine  the  natural  resonances 
and  filter  characteristics  of  the  entire  propulsion  system 

(2)  CEPSTRUM-like  analysis  of  the  predicted  ship's  power  spec- 
trum outputs  to  find  relatively  invariant  harmonic  groupings 
useful  for  passive  sonar  classification  purposes 

3.  3 ACQUISITION  OF  MONITORED  DATA  AND  ENGINEERING  DATA  IN 

PREPARATION  FOR  MODELING  ACTUAL  SHIPS 

Preliminary  efforts  have  been  made  to  organize  the  assembly  and  col- 
lection of  monitored  ships'  sound  data.  A substantial  amount  of  instrumenta- 
tion quality  recordings  of  ships'  radiated  sound  taken  under  controlled  condi- 
tions will  be  required  to  verify  the  computer  simulations  generated  through 
the  modeling  technique.  All  of  the  inherent  problems  of  acquiring  good  under- 
water acoustic  data  prevail  here. 

At  the  same  time,  more  than  one  promising  source  of  good  data  exists 
and  has  been  made  available  to  Litton.  There  remains  primarily  the  problem 
of  exploiting  these  sources,  which  will  require  that  time  and  a continuing  at- 
tention to  detail  be  devoted  to  this  aspect  of  the  work. 

As  a result  of  ship  monitoring  runs  carried  out  by  the  Navy's  Acoustic 
Pressure  Check  Range  (APCR),  Fort  Story,  Virginia,  tape  recorded  data  are 


in  hand  for  three  U.  S. 

Navy  ships: 

a. 

USS  Raleigh  (LPD  1) 

b. 

USS  Carp 

(SS  338) 

c. 

USS  Eaton 

(DD  510) 

Data  for  nine  merchant  types  were  gathered,  namely: 

a.  MS  Vishva  Shanti  (oil  engine) 

b.  SS  Drepanon  (reciprocating  steam) 

c.  SS  Queensville  (geared  turbines) 

d.  Copiapo  (geared  turbines) 
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e.  MS  Kerkedyk  (oil  engine) 

f.  SS  Shell  Naiguata  (geared  turbines) 

g.  MS  Hickory  Knoll  (oil  engine) 

h.  SS  Cosmopolitan  (turbo-electric) 

i.  MS  Athelbeach  (oil  engine) 

Three  of  the  ships  listed  above  were  visited  for  purposes  of  obtaining 
detailed  engineering  data.  These  were  Queensville,  Cosmopolitan  and  Athel- 
beach. Litton  observers  were  present  while  these  ships  were  being  monitored 
and  the  Litton  observers  considered  the  runs  to  be  sufficiently  free  of  inter- 
ference to  be  potentially  useful.  Since  then,  however,  the  taped  data  have  been 
submitted  to  STIC  for  spectrum  analysis.  It  appears  that  the  tape  recorder  at 
APCR  introduced  such  intense  artifacts  at  60  cps  and  harmonics  of  it  that  the 
data  may  not  be  useful.  However,  this  remains  to  be  investigated. 

Because  the  next  step  in  developing  the  ship  vibration  modeling/ sound 
field  synthesis  capability  is  to  apply  it  to  an  actual  ship  on  which  recorded 
sound  data  is  available,  the  SS  Athelbeach  will  be  modeled.  Excellent  engi- 
neering data  were  obtained  from  the  visit  to  the  ship,  and  additional  data  which 
will  support  a more  accurate  model  are  being  requested  of  the  ship's  builder, 
Hawthorne  Leslie,  Ltd.,  of  Newcastle  on  Tyne,  England.  The  ship's  engineer- 
ing plant  consists  of  an  oil  engine  (slow  speed,  two-cycle  diesel)  and  direct 
drive.  The  model  can  be  set  up  using  straightforward  modifications  of  the 
basic  model  which  is  presently  being  exercised. 
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SECTION  IV 

PLANS  FOR  FUTURE  EFFORT 


During  the  next  quarter,  effort  will  be  directed  toward  processing  re- 
corded ship  sounds,  continued  testing  of  processing  methods  on  simulated  sig- 
nals, an  analytical  study  of  the  important  problem  of  harmonic  group  detection, 
collection  of  additional  recordings  and  data  as  required,  and  a continuation  of 
modeling  the  mechanisms  of  ship  sound  generation. 

In  particular,  it  will  be  possible  to  determine  the  need  for  additional 
data  and  sound  recordings  only  after  analysis  of  the  initially  collected  record- 
ings. For  example,  it  is  not  clear  that  the  recordings  collected  at  Fort  Story 
are  adequate  for  a realistic  study  of  signal  processing  problems.  Further- 
more, it  is  desirable  that  several  examples  of  signals  of  significant  military 
interest  be  examined. 


4-1 


UNCLASSIFIED 

CONFIDENTS: 


UNCLASSIFIED 


REFERENCES 


(1)  J.  W.  Cooley  and  J.  W.  Tukey.  "An  Algorithm  for  the  Machine  Calcula- 
tion of  Complex  Fourier  Series,  " Mathematics  of  Computation,  Vol.  19, 
pp.  297-301  (April  1965). 

(2)  Christopher  Bingham,  M.  D.  Godfrey,  and  J.  W.  Tukey.  "Modern  Tech- 
niques of  Power  Spectrum  Estimation,  " Unpublished.  (Handout  and  talk 
at  the  meeting  of  Acoustical  Society  of  America,  Boston,  Spring  1966.  ) 

(3)  W.  M.  Gentleman  and  G.  Sande.  "Fast  Fourier  Transforms  — For  Fun 
and  Profit,  " Proceedings  of  the  Fall  Joint  Computer  Conference,  1966. 

(4)  T.  G.  Stockham,  Jr.  "High-Speed  Convolution  and  Correlation,  " Pro- 
ceedings of  the  Spring  Joint  Computer  Conference,  1966. 

(5)  R.B.  Blackman  and  J.  W.  Tukey.  The  Measurement  of  Power  Spectra, 
Dover,  New  York,  1958. 

(6)  W.B.  Davenport,  Jr.,  and  W.  L.  Root,  Random  Signals  and  Noise, 
McGraw-Hill,  New  York,  1958,  pp.  107-108. 

(7)  J.  P.  Den  Hartog.  Mechanical  Vibrations  — Fourth  Edition,  McGraw- 
Hill,  New  York,  1956,  p.  176. 

(8)  T.  P.  O'Brien.  The  Design  of  Marine  Screw  Propellers,  Hutchinson 
and  Company  Limited,  London,  1 9 62^ 

(9)  S.  Archer.  Torsional  Vibration  Damping  Coefficients  for  Marine 
Propellers,  Engineering,  1955. 

(10)  J.  P.  Den  Hartog.  Mechanical  Vibrations  — Fourth  Edition,  McGraw- 
Hill,  New  York,  1956,  p.  209. 


R - 1 


UNCLASSIFIED 


